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We study the contribution to vacuum decay in field theory due to the interaction between the 
long and short-wavelength modes of the field. The field model considered consists of a scalar field of 
mass M with a cubic term in the potential. The dynamics of the long-wavelength modes becomes 
diffusive in this interaction. The diffusive behaviour is described by the reduced Wigner function 
^ ■ that characterizes the state of the long-wavelength modes. This function is obtained from the whole 

I Wigner function by integration of the degrees of freedom of the short-wavelength modes. The 

dynamical equation for the reduced Wigner function becomes a kind of Fokker-Planck equation 
which is solved with suitable boundary conditions enforcing an initial metastable vacuum state 
trapped in the potential well. As a result a finite activation rate is found, even at zero temperature, 
for the formation of true vacuum bubbles of size M _1 . This effect makes a substantial contribution 
to the total decay rate. 
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I. INTRODUCTION 



In this paper we report our preliminary findings within a larger program which aims at the development of a 
theory of non equilibrium first order phase transitions, such as have occurred in the Early Universe (grand unified 
and electroweak symmetry breaking |jl|) and, possibly, in the first stages of heavy ion collisions (chiral symmetry 
breaking and confinement For this reason, we must seek a description of the decay process which emphasizes 

the dynamical aspects, over the static aspects encoded in the effective potential. 

Vacuum decay in field theory is described with a potential which displays a local minimum, separated from the 
absolute minimum by a potential barrier. A system prepared in the false vacuum state (metastable phase) within the 
Qh. potential well may decay in essentially two different ways, namely (a) by tunneling effect, that is, by going through the 
barrier in a classically forbidden trajectory [||-|(|, or else, (b) by activation, that is, by jumping above the barrier 
In either case, the decay probability follows the law P ~ Aexp(-B) which gives the probability per unit time and 
unit volume to nucleate a region of the stable phase within the metastable phase. In the tunneling effect, B = S/Ti, 
where Ti is Planck's constant and S is the action for the trajectory which goes under the barrier in imaginary time 
||. In thermal activation, B = V s /{kBT), where ks is Boltzman's constant, T the temperature, and V s is the height 
of the free energy measured from the false vacuum [^0|-p^[. Thus, activation disappears as T — > 0. 

In systems with few degrees of freedom, there must be an external agent, typically a thermal source, for activation 
to be possible. Our thesis is that in field theories there is a phenomenon similar to activation even at zero temperature. 
This comes from the observation that in a field theory, when a mode decomposition of the field is made, there are only 
a few long wavelenghts modes which are unstable and decay. These nearly homogeneous modes may be regarded as 
an open system, which interacts with the environment provided by the shorter wavelength modes. It is then possible 
to describe the quantum evolution of the system in terms of an effective dynamics, whereby the interaction with the 
environment results in the onset of dissipation and noise. 

The ultimate reason for the presence of a finite activation rate even at zero temperature is that, for a generic 
field theory the dynamics of these homogeneous modes is anharmonic enough to contain Fourier components with 
frequencies above the threshold for excitation of the short wavelength modes. This results in an energy transfer 
from the long wavelength or homogeneous modes to the short wavelength or inhomogeneous modes through particle 
creation. As demanded by the energy balance, and encoded in the fluctuation-dissipation theorem, this energy flow 
is compensated by a stochastic force on the homogeneous modes, originated in fluctuations of the inhomogeneous 
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modes. Thus the dynamics of the homogeneous modes becomes diffussive, even if, properly speaking, there is no 
external "environment" to the field |l3| . 

We wish to stress that this is not only a theoretical possibility. In this paper, we will show through a detailed 
analysis of a concrete model that the activation rate makes a substantial contribution to the full decay rate even at 
zero temperature. In the process, we shall develop the necessary formalism to compute the activation rate to leading 
order in h. 

The key ingredient will be the description of the quantum state of the long wavelength modes of the field by means 
of the reduced Wigner function. This function has the same information that the reduced density matrix of an open 
quantum system but is similar in many aspects to a distribution function in phase-space. The dynamical equation 
for the reduced Wigner function (master equation) includes noise terms produced by the short wavelength modes to 
quadratic order in the parameter coupling the short and long wavelength modes. We shall derive the tunneling rate 
from an analysis of the decay of non trivial solutions of the master equation, after suitable boundary conditions have 
been enforced. 

The master equation contains all required information on vacuum decay, including both the "tunneling" aspect of 
the homogeneous mode as well as "activation", i.e. the effect due to the back reaction of the inhomogeneous modes, 
but in this paper we concentrate in the activation aspect. Thus, whereas the instanton method ||, for instance, seems 
to be best suited to describe the pure tunneling of the homogeneous mode, the effect due to the back reaction of 
the inhomogeneous modes seems to be best described by the present method. Another advantage of this approach, 
which was one of our first motivatons, is that whereas the instanton method works well for processes not far from 
equilibrium, this method should work also for out of equilibrium situations since it is based on a real time evolution 
equation such as the master equation. 

Our conclusion shall be that, for a physically well motivated system-environment separation, the contribution to the 
exponent coming from quantum fluctuations in the short wavelength modes is comparable to, but distinguishable from, 
the contribution from tunneling under the barrier by the long wavelength modes themselves. Which contribution is 
actually dominant will depend on the specifics of each model. Therefore, the approaches that underplay backreaction 
from the short wavelenght modes, underestimate the tunneling rate, and may even miss the largest contribution. 

Comparing with the literature on tunneling in open and driven systems, the main difference is that we do not assume 
a priori any features of the noise and dissipation, but rather derive them from the underlying unitary field theory. 
In practice, this means that we shall have to deal with non-local dissipation and coloured noise. Also for simplicity, 
we shall focus on computing the exponent in the tunneling rate to leading order in ?i, and at zero temperature. As 
a matter of fact, this is the difficult case, since at high temperature the interaction of the short and long wavelength 
modes is just another aspect of the interaction with the heat bath. 

The paper is organized as follows. In Sec. |n] we briefly review current theories of vacuum decay and place our 



work in that context. In Sec. Ill we present the basic derivation of the master equation for the reduced Wigner 



function for a Quantum Brownian Motion (QBM) model which is of relevance for our problem. The model is an open 
quantum system consisting of of a massive particle in an arbitrary potential coupled to an infinite set of harmonic 
oscillators, the coupling is linear in the system variables but quadratic in the oscillator coordinates. For our problem 
the master equation reduces to a Fokker-Planck like equation. In Sec. [w] we consider a field theory model consisting 
of a massive scalar field with a negative cubic potential term. The long and short wavelength modes of the field are 
separated and we reduce the problem to a system-environment interaction similar to the QBM model discussed in the 
previous section. The relevant kernels for the corresponding Fokker-Planck equation are calculated. Sec. |v| contains 
the derivation of the tunneling rate from the analysis of the solutions to the Fokker-Planck equation. We discuss our 
results in the final section. A number of appendices carry some of the unavoidable technical details. 



II. VACUUM DECAY: A SURVEY 



1. Vacuum decay in systems with one degree of freedom 

The twin issues of thermal activation and spontaneous nucleation have a long and distinguished history; see Refs. 
Jl4] , [ro| for a review. In its simplest formulation, we deal with a quantum mechanical system with one degree of freedom 
x and Hamiltonian H = p 2 /2 + V{x), where V has the generic form shown in the figure. It is an experimental fact 
that even if we prepare the system to be confined to the neighborhood of the "false vacuum" x ~ 0, and the barrier 
is much higher than the typical energies accesible to the system (E ~ KV" (0) in the quantum mechanical problem, 
E ~ kT at finite temperature), the system will find a way to escape from the potential well after some typical time r 
has ellapsed. The problem is to estimate the "mean life" r or equivalently the "tunneling rate" r _1 . 
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Vacuum decay can be formulated as a steady state problem if we inject particles into the system in order to keep 
a constant population in the false vacuum state. We then have a constant flux of particles impinging on the barrier 
from the left, and the problem reduces to the computation of the transmission coefficient In quantum mechanics, 
this is readily obtained within the WKB approximation, and the result is the Arrhenius-like expression 

t- 1 ~ Acxp(-5/fi), S= dx y/2V(x), (I) 

JO 

where A is a prefactor of order 1, and we have set the classical energy of the false vacuum to 0. 

This formula describes vacuum decay through tunneling, that is, an essentially quantum phenomenon. If we allow 
the system to interact with an external noise source (typically a heat bath at a given temperature T), then the energy 
of the system alone is no longer conserved, and the system can jump over the barrier, resulting in vacuum decay 
through activation. The activation rate has been computed, within the "constant flux" approach, by Kramers m and 
Langer M. They show that equation (|l|) still holds, but the exponent S/U becomes F s /kT, where F s is the activation 
free energy, i.e. the height of the free energy barrier to be overcome. 

A less artificial approach attempts to compute the actual time evolution of the false vacuum state f p. This is 
not a stationary state of the system, but it may be expanded in energy eigenstates tpE- The amplitude of i^e hi the 
expansion of ^fp peaks around E ~ 0, and for small energy may be approximated by the Breit-Wigner form [^[[j) 

ff,(*,t)~ JdEe^p^. (2) 

At long, but not too long, times we obtain * F - Ce- xt / n . The false vacuum behaves as an energy eigenstate 
with complex energy Ep ~ —iX (complex energies and eigenstates may be defined within an extended formalism of 
quantum mechanics Jig] ), and the mean life is r ~ HX^ 1 . 

The calculation of activation rates in the "complex energy" approach has been worked out by Langer Jig] . The 
idea is to define a "free energy" for an ensemble of unstable configurations including the critical droplet. This "free 
energy" is complex, and its imaginary part is related to the mean life as in the quantum mechanical problem. The 
physical free energy, of course, is real, and it is given by Maxwell's construction pj^ . 

Coleman and collaborators have proposed a simple and elegant way to compute the complex false vacuum energy 
. The idea is that the vacuum energy can be expressed in terms of a path integral over Euclidean histories with 
appropiate boundary conditions. For unstable systems, the path integral must be computed by analytic continuation, 
and an imaginary part appears. In certain cases it is possible to show that the path integral is dominated by the 
contributions from a discrete set of saddle points, corresponding to sequences of "bounces" against the inner sides 
of the barrier with little or no overlap between bounces. Then the formula in Eq. (Q) is recovered, where S is now 
interpreted as the euclidean action for the one bounce solution, also called the "instanton" . 

The instanton method is easily generalized to the thermal case jll] . The idea is to write the partition function for 
the unstable system as a path integral over Euclidean configurations with periodicity (5% (/3 = 1/kT) in Euclidean 
time, and then to evaluate the path integral in the saddle point approximation. Due to new boundary conditions, the 
thermal instanton may not be the same as the instanton at T = 0. The change in the nature of the instanton gives a 
simple and compelling interpretation of the crossover from spontaneous transition to thermal activation p"2"[ . 
The tunneling rate can also be derived from the large order behavior of perturbation theory [ fl8| 



2. Tunneling in systems with few degrees of freedom 

All approaches discussed so far have natural generalizations to systems with few degrees of freedom. In the case of 
the instanton approach, the generalization is almost immediate, only one has to take care of symmetries of the system 
which may appear as zero modes in the spectrum of perturbations around the instanton solution, thus causing an 
apparent divergence of the path integral. These symmetries may be handled by isolating them as collective modes 
prior to the saddle point evaluation of the path integral . 

The "constant flux" approach is implemented by seeking a solution to the Schrodinger equation within the WKB 
or Born-Oppenheimer approximation pc|[ . The idea is to identify a single variable x which parametrizes the "most 
probable scape path" , namely the path across the saddle separating the false and true vacua. Then one uses a mixed 
ansatz for the wave function, whereby it is assumed to be Gaussian on all other variables, and of WKB form with 
respect to x. "Under the barrier" , the WKB approximation leads to a Hamilton-Jacobi equation for a particle moving 
in an inverted potential, the same dynamical problem one confronts in the instanton method, although now the 
potential may be modified by terms arising from the zero point energies of the transversal modes; but not of x itself 
pl[. This solution is then matched to oscillatory solutions on either side of the barrier. 
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Conceptually, the WKB method has some advantages over the instanton method p2[ . In the first place, the 
connection to the Schrodinger equation is much more straightforward. Also, one has more control on the quantum 
state of the transversal degrees of freedom, which then allows one to ask questions like whether tunneling is associated 
to particle creation p3| - ^of. Finally, it has proved easier to define the range of validity of the approximations involved 
in the WKB framework than in the Euclidean path integral one. However, they are likely to be equivalent in last 
analysis, and the conclusion that these methods work best when the tunneling variable x is slow compared with the 
transversal modes probably holds equally well for both approaches |2f|. We refer the reader to Ref. |27j for recents 
developments along these lines, and to Ref. [p8| for an example of similar techniques in a different context. 

Some nonequilibrium aspects of vacuum decay have been the focus of work by Boyanovsky et al. j2j|. They show 
that realistic initial conditions usually imply a nonzero probability for the system to be at the unstable configuration 
in the saddle of the free energy surface. Starting from this configuration, the roll down of the system towards the 
true vacuum may be analyzed by usual nonequilibrium field theory methods. However, the initial amplitude is still 
computed by conventional methods, such as discussed so far. 



3. Tunneling in quantum field theory 

Gcrvais, Sakita and De Vega have applied the WKB method to tunneling in quantum field theory see Refs. 
|3l| , j32f for instanton and complex-time methods. In spite of the obvious similitudes, there are some important 
differences between the problem of tunneling in field theory and in systems with few degrees of freedom. Some of 
these differences are technical in nature, such as the need to carefully account for loop corrections to the effective 
potential |33[ , and to adopt a regularization procedure to compute the prefactor in the Arrhenius formula [^-^6), 
which would be divergent if naively computed. There is also a fundamental conceptual difference which we now 
discuss. 

A field theory only makes sense, from the physical point of view, if it is understood as an effective theory describing 
the low frequency sector of a more fundamental theory, whose high frequency degrees of freedom may be totally unlike 
continuous fields; such as strings in elementary particle physics, discrete lattices in condensed matter applications, 
molecules in hydrodynamics [p7[ . 

One clear way to bring this point home is by explicitly integrating out all modes with wave number k > A, where 
A is some cutoff. The remaining modes are described by a coarse grained effective potential p8|-[43f . As A decreases 
from oo to 0, the coarse grained effective potential interpolates between the "bare" potential shown in the figure to 
the Maxwell construction free energy , showing the drastic effect of the short wavelength modes or high frequency 
sector on the physics of the long wavelength modes or low frequency sector. At some point the barrier separating the 
different metastable points dissapears, reflecting the effect of averaging the field over distances much larger than a 
domain. 



4- Tunneling in open systems 

As we have seen from earlier studies of the coarse grained effective potential, tunneling in field theory should be 
properly posed as an open quantum system problem. However, the nature of tunneling in an open system cannot 
be described adequately by just computing changes in an effective potential. Besides the static changes reflected by 
the scale dependent effective potential, the dynamics of the long wavelength modes will become both dissipative and 
stochastic. The onset of dissipation and noise is also generic to semiclassical [|l3j and effective theories j4f|. 

For the present discussion it is essential to realize that noise and dissipation are actually two aspects of a single 
phenomenon, the dynamical action and back-reaction between "system" and "environment" . In equilibrium situations, 
this inner relationship can be made explicit through the "fluctuation-dissipation" relation JfC] ] . 

A simple way to deal with tunneling in open systems is to model the environment explicitly within a larger system- 
environment complex, in effect reducing the problem to tunneling in many degrees of freedom p7|-|49|]. However, it is 
essential to avoid approximations (such as assuming that the environment degrees of freedom perform linear oscillations 
around a prescribed trajectory of the system) which in practice underplay the back reaction of the environment on the 
system, and thus break the balance between fluctuation and dissipation. The relevance of the fluctuation-dissipation 
relation to tunneling has been emphasized in refs. 15(3]. 

The broadening of the reduced Wigner function of the open system by external noise has been discussed in Ref. 
pl| . Unlike the present work, these authors consider an external noise source, whose spectral features may be chosen 
at will. In order to give an adequate account of back reaction, Bruinsma and Bak |p2| have proposed treating the 
system as propagating in a random medium, the randomness being associated to the environmental variables. In a 
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second step, a path integration over histories of the bath allows the computation of the tunneling rate. As in the 
present work, a serious consideration of back reaction leads to describing the system as a driven system, subject to 
stochastic forces originating from the environment. 

The theory of vacuum decay in open systems has points of contact with the problem of decay in driven systems 
p3| , E4j , although in these later studies usually the properties of the driving force are assumed a priori, rather than 
derived from a more comprehensive model. It is also possible to obtain a path integral representation of the solution 
of a Langevin equation, whereby an open system may be subsumed into a larger field theory [ fi5| , [56| . See Ref . |5^] for 
further developments. 



5. Tunneling and the semiclassical approximation 

Semiclassical field theories, where some fields are treated as c-number, while the rest is described quantum mechan- 
ically, may be seen as open systems, with the classical fields as the system and the quantum fields as the environment. 
Quantum fluctuations in the system are registered as noise by the environment, and may induce transitions. 

An early application of these ideas appeared within Starobinsky's "stochastic inflation" program (58). Here the 
superhorizon modes of the inflaton field during inflation are the system whereas all other shorter wavelength modes 
are the environment; for a discussion of the validity of the semiclassical approximation, see [ p9| , p0[ . Cosmological 
redshift causes a continuous streaming of modes from the environment to the system, which may be regarded as a 
white noise source. This noise may allow the system to hop over potential barriers, seeking the absolute minima. 
Eventually, a steady state distribution of cosmological domains is reached, not unlike that predicted by the Hartle 
and Hawking "wave function of the Universe" ]22|,[3l| . 

Much more generally, within the semiclassical approximation the back reaction to cosmological particle creation 
processes always has a stochastic component [[l3| , for which reason the correct semiclassical description of the Early 
Universe ought to be formulated in terms of an stochastic "Einstein-Langevin" equation ]62| . If we consider an 
ensemble of Universes, then we may introduce a distribution function obeying a Fokker-Planck-like equation |]63|] . This 
equation describes activation phenomena, which are the semiclassical version of Vilenkin's "creation from nothing" 
scenario Q. 



III. OPEN SYSTEMS AND THE REDUCED WIGNER FUNCTION 



Before dealing with field theory we will consider in this section a Quantum Brownian Motion (QBM) model which 
is typically used as a a paradigm of an open quantum system. The system has an arbitrary potential and the coupling 
between the system and the environment is linear in the system variables but quadratic in the environment variables. 
That feature will be of relevance when dealing in section IV with our field-theory model. The main result of this 
section is the derivation of the master equation for the reduced Wigner function of the QBM model to leading order 
in H. This equation turns out to be a Fokker-Planck equation which is similar to that used by Kramers |7j to study 
the activation problem in statistical physics. 



A. A QBM model 



As our QBM model we consider a system consisting of a particle of unit mass (M = 1) described with a variable 
x and subject to an arbitrary potential with a quadratic part corresponding to an oscillator of frequency f2o and an 
anharmonic part V^ nl '(x), i.e. V(x) — {1/2)0^ x 2 + V^ nl '(x), which is coupled to an environment consisting of an 
infinite set of harmonic oscillators with coordinates qj . The action for the whole set of degrees of freedom is defined 
by 

S[x,{qj}] = S[x] +S[{qj}} + S* mt [2, {<?.,■}], (3) 
where the system, environment and interaction actions are given, respectively, by 

S[x] = fdt (\x 2 - V(x)) , (4) 



(•5) 
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S lnt [x, {q 3 }} = ^g dt xq?, (6) 
3 J 

where g is a coupling constant, and we have assumed that the coupling is linear in the system variable but quadratic 
in the environment variables. The environment oscillators have all the same mass m and their frequencies are ujj. At 
this point the potential V ljll \x) is arbitrary, but later we will take a cubic potential, V^ nl '(x) = — (A/6)x 3 , in this way 
the total potential will present a local minimum and a barrier as required to represent the system in a metaestablc 
phase. Also the parameters A and g are unrelated, however when we consider the application to a field theory these 
parameters will coincide. 

The reduced density matrix for our open quantum system at a certain final time tf is defined from the density 
matrix p of the whole system by tracing out the environment degress of freedom at that time 

p r (x f ,x' f ,tf) = l'[\_dqjp{xf,{q j },x' } ,{q j } i tf). (7) 
J 3 

The reduced density matrix at time t / can be written in terms of the reduced density matrix at the initial time ti 
by the evolution equation, 



p r (xf,x'f,tf) = J dxidx' i J(xf,x' f ,tf,Xi,x' i ,t i )p r (x i ,x' i ,t i ), (8) 
in terms of the propagator J, whose path integral representation is 

J(xf,Xf,tf,Xi,x i ,ti)= / Vx / Vx' exp — (S[x\ — S[x'] + Sif[x, x']) . (9) 

J Xi Jx[ 

where Xi = x(U), Xf — x(tf) and similarly for the primed quantities, and Sif[x, x'] is the Feynman and Vernon 
influence action [ p5[ . When the system and the environment are initially uncorrelated the initial density matrix 
factorizes, i.e. p(Ji) = Pr{U)Pe(ti) (here p e stands for the environment density matrix), the influence functional, 
which is defined by F[x,x'] = exp(iSiF[x, x']/h), can be expressed by 

r rq IJ) rq u) 

F[X, X'] = dq^dqfdq'V J ' Vq 3 j ' Vq' 3 exp - (S[{q }\ - S[{q' 3 }\ + S int [x, {%}} - S int [x' , {g}}]) 

Mi^hWPhU), (10) 

where q^ — qj(U), q'^ = q'jiti), and at the final times qj(tf) = q^ = q'j(tf). 

Assuming a Gaussian initial state for the environment, p e is Gaussian and the influence functional can be computed 
perturbatively in g from the path integral. Up to second order in g [p{| we have for the influence action, 

*/ 



Si F [x,x'] = -2 / dt j dt'A(t)D(t,t')X(t') + - j dt j dt'A(t)N(t, t')A(t'), (11) 



where we have introduced the average and difference coordinates defined, respectively, by 

X(t)=±(x'(t)+x(t)), A(t)=x'(t)-x(t). (12) 

The kernels D(t,t') and N(t,t') are called dissipation and noise kernels, respectively, and are defined by D — J2j Dj 
and TV = J2j Nj where 

D j (t,t') = - l -{[E j (t),E j (t')}), (13) 



Njfrt) = -({Ejit),^')}) - (E^iE^t')), (14) 



with Ej = gqj. It is now convenient to introduce the kernels Hj(t,t') = —2Dj(t,t')6(t — t') and we can write the 
influence action in the form 



G 



Si F [x,x'] = A • H ■ X + |A • N • A, 



(15) 



where we have introduced the notation A ■ B = J dtA(t)B(t) and defined H — Hj which we may write formally 
as H(t,t') = —2D(t,t')6(t — t'). This last equality is however a formal expression since being the product of two 
distributions, H is not well defined and suitable regularization and renormalization is required. This term, in fact, 
has local divergent parts that may be reabsorbed into the parameters of the bare action, see p6| for details. Thus, 
from now on we will assume that if is a well defined distribution in the previous sense. 



B. The reduced Wigner function 

Our main purpose in this subsection is to write the reduced Wigner function for the system in a suitable way. This 
is a phase space function which is defined from the reduced density matrix by the following integral transform 

f{X ^ t)= 2^hj dA e ipA/h Pr (X-A/2,X + A/2,t). (16) 

The reduced density matrix (||) at time tf can be computed from the path integrals of (||). To carry out this 
computation we will follow closely Ref. J67| where a similar computation for a linear system coupled linearly to an 
environment was described. For this reason we will describe here only the main steps and will concentrate on those 
which are peculiar to the nonlinear system. This is performed in several steps in which a key role is played by the use 
of the coordinates X and A instead of x and x' . The first step is to integrate the system action in (^) by parts using 
the new coordinates X and A, 



ftf / j2 pfw 

S[x] - S[x'] = -XA\l( + J dt A(t) (^X(t) + M 



A=0 



(17) 



where the ellipsis stands for the terms nonlinear in A that come from the potential V{x) which involve higher 
derivatives in X (see below). Note that due to the fact that the potential gradient is evaluated at A = this term 
can also be written as V'(X) where V(X) is functionaly the same as ^(a;). 

The change of integration variables f* f Vx J*/ Vx' — > f* f VX J^ f T>A involves a Jacobian which is unity and 

thus the path integration of the propagator (j^) can be written as 

1 VX [ ' VAe^ A - L ^F[A,X], (18) 



where L[X] is a functional of X and a function of t defined by 

dV x rt 



L[X;t)^[- 2 X { t) + - 



A=0 



dt'H(t,t')X(t'), (19) 



and the functional F[A, X] incorporates in the exponent all the terms that are not linear in A which come from the 
influence action and from the nonlinear potential of the system action, when expressed in the variables X{t) and A(t). 
More precisely, 



F[A, X] = exp % - ( l -A ■ N ■ A + / dt' V (nl) [A, X] ) , (20) 



where 



V {nl) [A,X] = 2 



Q{2n+l)y 



^ (2n+ 1)! dX( 2n+1 ) 

n>l 



(-A/2) 2n+1 . (21) 

X,A=0 



In particular, for the cubic potential V^(x) = -(A/6)x 3 we have V^[A,X] = - (A/24) A 3 . 
Let us now introduce the functional Fourier transform 

P[i] = K (VA e <A.C/R e -A.JV.A/2R (22) 
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where K = 1/ det(27r7il), the interpretation of this functional will be discussed below. With expression ( p2| ) we may 
write the reduced Wigner function as 



VA e ,A.( i [i H )/ fie ./v«""/s e .A,x (i ,)/ Sfr(Ii _ At/2 ,X l + Ai/2,*,); 



A/ 
A, 



and using that A = exp(iA • £/K) (i%S/S£) exp(— iA ■ £,/h), it may be rewritten as 



f(X f)Pf ,t f ) =2^/ dXd^JJvX J V^P Q %X;t f ) J°° dA f e <Pf-*<.tf))&f/* 
VA e ^-{L[x]-i)/h e iA i x{n)jn p ^ x . A j 2 ,X t + A l /2,t t ), 



A/ 

A, 



where 



P Q [tX-t) = exp 



(23) 



A convenient way to perform the path integration for X(t) is to introduce the following functional change: 

X(t) - {x = ft = X(U),i(t) = L[X;t)} . (24) 

With this transformation the function X(t) becomes substituted by the initial conditions (Xi,pi) and the function 
l;(t) in the path integration. This functional change is invertible, in the sense that {Xi,pi,£(t)} — > X(t), since the 
solution X(t) of the integro-differential equation involved in ( |24| ) is unique given initial conditions (X^pi). A subtler 
point concerns the Jacobian of the transformation (|24|). Even though this transformation is nonlinear one can show 
that the Jacobian is constant. This can be seen by discretizing the time tk = ek + ti (k = 1, 2, . . . , n, and U is the 
initial time). Then the corresponding values (X,, Xi, . . . ,X n ), where Xk — X(tk), map into (Xj,p,, ■ ■ ■ , £n) m 
such a way that the Jacobian matrix has zero elements above the diagonal. For instance, the second derivative terms 
become [(Xk — Xk-i) — (Xk-i — Xk-2)]/^ 2 + V'(Xk-i) — Ylk'<k Hkk'Xk' = £fc. The Jacobian is thus the product of the 
diagonal elements, which are constant (independent of any Xk). Then one may write J VX . . , = K J dXidpi J Vt; . . . 
and introduce convenient delta functions such as S(X(tf) — Xf) to ensure that the correct final points appearing in 

(|l8|), i.e., J 1 VX are recovered from the functional integral J V£ with free ends. One should also be careful about 
the dependence on the initial conditions (Xi,pi) in the general case. 

Now we first perform the integral VA which simply leads to a term proportional to <5(£ — £) and the integral 2?£ is 
then trivial. On the other hand the integral dAi brings back the reduced Wigner function at the initial time according 
to Eq. (0). Finally, we get the following suggestive form for the reduced Wigner function at the final time 

f(Xf,p f ,t f )=K J dX t d Pl J V£ P Q [£;t f ) 6 (p f - X (t f )) 6 (X f - X (t f )) fiX^p^U), (25) 

where X(t) is a solution of the integro-differential equation 

L[X;t) = ft), (26) 

with initial conditions (Xi,pi), i.e. X = X[£; Xi,pi), and 

P Q [t;X i ,p i ,t)=P Q [t,X[t;X i ,p i );t). (27) 

The constant K form the Jacobian can be determined from the condition that Trp r (tf) — 1, i.e. that 
J_ dXfp r (Xj, Xf, tf) — 1, which is equivalent to dXfdpff(Xf,pj,tf) = 1. Inserting expression (25) for 
f(Xf,p f ,tf) we get 

K J dX. t dp t J VPQfctfWXi^U) = 1. (28) 
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When Pq[£,X;t) does not depend on X, one can use the fact that dXi dpj/(JQ,pj, i,) = 1, which is as a 
consequence of Tr p r {ti) = 1 to obtain, 

K Jv^P Q [^;t)=K Jv^P[(\ = l, (29) 

which determines K. 

Several remarks are in order here. The functional Pq[£, X[£; Xi,pi)\t) is always real, but in general it will not be 
positive definite and, thus, will not really correspond to the probability density functional for a classical stochastic 
process. This is the meaning that must be associated to the stochastic process in the Langevin-like equation (p(j|). 
This situation is, in fact, analogous to that for the Wigner function but applied here to distribution functionals. 
Note that this is in contrast to the linear case studied in Ref. j6?|] , where the source of the Langevin equation really 
corresponded to a stochastic process (with a positive probability density functional). 

We emphasize again that if we have a cubic potential for the system, V(x) = — (A/6)cc 3 , and keep up to quadratic 
order in g, we have explicitly 

F[A,X]=e-^ A - N - A e~U dt ^ A \ (30) 
where the noise kernel N and the kernel H, which appeared in L above, are both quadratic in g; note that there is 



no dependence on X in this case, and thus Pq[£;£) defined in ( |27[ ) does not depend on the initial conditions (Xi,pi). 
Hu, Paz and Zhang ]68[ | obtained the master equation for the particular case in which the nonlinear potential of the 
system is also treated perturbatively in A, which was considered to be of the same order as g. Here, however, our 
result is exact in A, and to leading order in fi. This fact will turn out to be important since in Sec. |fv| it will be 
crucial to consider solutions of the classical equations of motion which are nonperturbative in A thus reflecting their 
strong nonharmonicity. 



C. The master equation 



The expression (p9) of the reduced Wigner function and the Langevin-like equation ( pq ) can be used to derive the 
master equation for f as a formal Fokker-Planck equation. The derivation is usually handled using Novikov's formula 
when the stochastic process is Gaussian |32],[73 67 1. Here, however, this is not the case for Pq and we have to work 
from the beginning. 

To obtain the equation of motion for the Wigner function, we derive both terms of eq. ( p5| ) with respect to time. 
Observe that Pq depends explicitly on time, therefore 



-^f(X f ,p f ,t f )=A + B, 



where 



(we write t f in Pq to emphasize that the dependence on t is taken care of explicitly by the B term) and 



A = K dX,d Vl / VSP Q [£;fj) 



p f -X(t f ))S(X f -X(t f )) )f(X i , Pi ,t i ) 



B = K 



Let us analyze the B term first. Since 



integrating by parts we find, 

B = K f dX,d Pl [ V£P Q [S;t f ) 



k y(nl) [lh satir Xf] ( 5 ( p/ " * 6 (Xf ~ x {tf}) 



(31) 



(32) 



J dXidpi J Vi (JLpyfctffj d{p f -X{t f ))5{X f -X{tf))f{X u p h U). (33) 



(34) 



f(X llPl ,U) (35) 
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We are only interested in derivatives taken at tf, when 



SX(tf) 



0, 



= 1, 



(36) 



without further dependence on This can be seen from the fact that 5X(t)/6£(t') satisfies ((L + H) ■ 

SX I '6£)(t,t') — S(t — t') with the X(i) which appears in L fixed, L is the integro-differential operator L(t,t') = 
(d 2 /dt 2 + d 2 V/dX 2 \ xm ) S(t - t'); see Eq. ©. The solution is 5X(t)/S^(t') = G ret (t,t'), which is the retarded 
propagator corresponding to the linear operator (L + H) with X(t) fixed. So the final result is 



B = -vM[-ih—,X f ]f(X f , Pf ,t f ). 



Concerning A, we find 
A = -K f dXidpi [ V£ Pg[t;tJ\ 



dX(t f ) d dX{t f ) d 



8[pf-X{tf))8(X f -X(tf)) 



dt f dX f dt f dp f 

and reading the derivatives from Eq. ( |26| ) we can write A = A\ + A^ + A%. The first term is simply 

Ax = {H s J}, 



(37) 



(38) 



(39) 



where {H s ,f} = —p(df/dX) + V {df /dp) is the Poisson bracket with H s = p 2 /2 + V(X) the system Hamiltonian. 
The second term is 



d 
dp/ 



K / dX l dp l / VtP Q [btl) 



dt'H(t f ,t')X(t') 



s(p f -X(t f ))s(X f ~X(t f )) 



fi ) 



(40) 



to lowest order in h, we are entitled to replace X(t') inside the non-local term by a solution of the classical equations of 
motion with the given Cauchy data Xf and pj. We shall call this procedure of substitution of the classical trajectories 
into the terms which are already of order h "reduction of order" . We may then extract the non-local term from the 
integral to get the simpler form 



a 2 = ^-[r(x,p,t)f], 

where T{X,p,t) = -J* dt' H(t,t')X(t'). Finally, the third term is 



A 3 - - Ajr j dX t d Pl J Z>£ P Q [£ : tj) £ (t f ) S(p f -X (t/)) 6(X f -X (t f )) f t . 



dp f 

To compute this term, we note from (^) and ( p7[ ) that 

£(tf)P Q [&tj) = L X p i f 

and since P[£] is Gaussian, we may use the identity (Novikov's formula): 

ti(tf)PQ[Z ■■ tj) = dt'N(t,t')-^j-P Q [^tJ) 



7 dt ' vW[-ih4z,X] 



(41) 



(42) 



(43) 



(44) 



Integrating Eq. ( ^2|) by parts and after further simplification, where we explicitly assume that, as in the case of 
the cubic potential, Pq is independent of X (in this way we can commute the exponential of in Pq with the 

functional derivative with respect to £), we obtain 



A3 = ^{NJ}, 
dp f 



(45) 



where N = hj dt' N (t,if) X (t 1 ), and we have applied once again a reduction of order prescription. The details of 
this calculation are given in Appendix A. 
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To summarize, and using the explicit form of for the cubic potential, we obtain the following dynamical 

equation for the reduced Wigner function (master equation): 

% {"-/I ' i:im,P,t)f} + ^-{NJ}~h 2 ^-^r. (46) 



dt ' dp dpf ' 24 dp 



If the system were isolated, the master Eq. (46) would reduce to 



dW o A d 3 W 

_ = W ,^_ S *__ (17) 

where W is the Wigner function of the closed system, H s its Hamiltonian and the curly brackets are the Poisson 
brackets. This equation is exactly equivalent to von Neumann's equation for the density matrix of a one-dimensional 
quantum mechanical system with a cubic potential V(x) — — (A/6)x 3 . Note that the term with the third derivative 
with respect to the momentum is the responsible for tunneling when properly combined with the otherwise classical 
dynamics generated by the term corresponding to the Poisson bracket, i.e., if this term were not present, the evolution 
of the Wigner function would be entirely equivalent to that of a classical ensemble in phase space. 

There is a theorem by Pawula |3J| which states that a diffusion-like equation such as Eq. ( |47| ) should have up to 
second order derivatives at most, or else an infinite Kramers-Moyal expansion, for non-negative solutions W(x,p, t) 
to exist. The equation for the Wigner function circumvents the implications of the theorem since it need not be 
everywhere-positive. Even if we have an everywhere-positive Gaussian Wigner function at the initial time, the 
evolution generated by an equation such as Eq. ( ^7|) will not keep it everywhere-positive. This can be connected 
with the fact that the Wigner function could be interpreted as the distribution function associated to an ensemble 
of solutions satisfying the Langevin-like equation ( p6|) with a generalized stochastic source £(£) with a non-positive 
probability density functional -Pq[£]- Thus, here we see the essential role played by the non-positivity of the Wigner 
function in a genuinely quantum aspect such as tunneling. In other aspects such as in quantum coherence this role is 



well known 70 



It is important to stress the following points. First, we are not assuming that the stochastic trajectories described 
by the Langevin equation are real trajectories; although they may be if there is decoherence ffl] ]. Second, although 
we use the initial Wigner function to weight the initial conditions, we do not asume it is a probability distribution 
function. Our interest in the Wigner function is that we shall use the fact that it acquires a substantial nonzero 
average beyond the basin of attraction of the false vacuum as a signal that tunneling has occurred. This application 
is valid even if the Wigner function itself cannot be understood as a probability distribution function, because the 
distributions computed from the Wigner functions, such as J^^dX f(X,p) or dp f(X,p) are true probability 
distributions, and dX (p/M)f(X,p) is a true probability flux |72]. 



D. Dynamics of the distribution function / 

To compute the tunneling probability from equation ( f47f ) is possible using, for instance, a WKB approximation 
scheme p0|j2~l| . Alternatively one may use the instanton method |)| which gives a simple answer in this case. Our 
main interest in this paper is not this contributon to tunneling but rather to compute the effect due to the back 
reaction of the environment. To simplify our derivation we shall assume that the time scales are different, that the 
dominant term is the back-reaction effect and, thus, we will neglect the third derivative term in the master equation 
([46]). More precisely, we will use the evolution equation for the distribution function / 

ft = {HsJ} + j^ x >^ + dT f {NJ} (48) 

This equation describes an ensemble of points evolving according to the dynamics 

X=p, p=-V'(X)+ f dt' H (t - t') X {t') + (49) 



where H is the dissipation kernel introduced in Eq. (|15|), and £ is a Gaussian noise described by the noise kernel iV 
introduced in (|Tl]), as 

(£(*)£(*')> = Mr (*-*')• (50) 
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Although we set up the initial conditions at ti = 0, we extend the lower integration limit in Eq. ( |49| ) to — oo 
for computational purposes. The approximation is, nevertheless, justified since the characteristic frequencies of the 
environment, u>j > 1 (recall that the system particle has unit mass), are much larger than the typical decay rate of 
the initial false vacuum state. In other words, the characteristic time scale for the environment dynamics t < 1 and 
hence, the typical "correlation time" for the dissipation and noise kernels, is much smaller than the typical decay 
time. The contribution to the integration interval (— oo, 0) is, therefore, relatively small. 



According to Eqs. (13) and (|14|) these kernels admit the following representation, 



N{t-t')=\ ({3 {t) , 3 (*')}) - (3 (t)) (3 (*')> , (51) 

H(t-t') = ~<[3(*), 3 (*')]> #(*-*'), (52) 

where 5 = y~]j Ej and where we must keep in mind the need for regularization of the kernel H as defined above. Let 
us take the Fourier transforms, 

<[3 (t) , 5 (*')]) = J ^ e-^-^h^ M , (53) 

N (t-t')= I — e-^-t'h (lu) . (54) 
J 2tt 

For an environment initially in thermal equilibrium at temperature T = /3 _1 the functions 7 and v will be related 
through the fluctuation-dissipation theorem 



l + fo H 



K M 7 H , f (u J ) = {e^-l) \ (55) 



which is a consequence of the KMS formula. Here we shall consider the zero temperature case only. 

In order to compute the memory terms in equation (^), it is convenient to parametrize the trajectories by their 
initial conditions at time U = 0. These trajectories may be written in terms of the action-angle variables J and 9 
associated to the classical potential, which we assume it has a potential well bounded by a finite potential barrier. 
In other words, we are using J and 9 as Lagrangian coordinates, identifying a given trajectory, while X and p are 
like Eulerian coordinates, identifying where the trajectory is at a given time. The action variable is defined by 
J = (l/27r) <j>pdX; since p can be written in terms of X and the system Hamiltonian H s , substitution into the 
equation defining J and inversion implies that H s — H S (J), and dH s /dJ — Q(J) is the frequency of the motion. The 
angle variable 9 changes from to 2-7T and satifies the equation of motion 9 = Q. 

Since the kernel H is already of order h, in the memory term of (^9] ) we must replace the trajectory X (t 1 ) by a 
solution of the classical equation of motion, in which case the transformation to variables (J, 9) is canonical. For fixed 
J, the classical trajectory is periodic, X(t') = X(t' + 2ir/£l); note that the motion is periodic inside the potential well 
but at higher energies, near the top of the potential barrier (the separatrix) when J — > J s , the motion ceases to be 
periodic and O — s- 0. Thus we may write 

X (t') = e m [ 9+0(J)t 1 X n (J) (56) 

n 

p{t') = 1 J](J)^e m [ e+!!(J)t ']jiX 11 (J) 

n 

where X- n = X*, since X is real. We then write the memory dependent term as: 

T = -f dt>H{t-t')X{t>) = -Y,X n {J)^ + ^ -gffij . , (57) 

J_oo J 27r oj + nil ( J) - ie 



where we have used (pj) , ( p3|) and (|56|), and that / °° dwexp(isu) = i/(s + ie). Therefore we have 

T(X,p,t) = -J2X n (J)e^Whn(J), 7n(J) = f P—^M— (58) 

Observe that although the Langevin equation is now local in time, it is not necessarily ohmic. A similar manipulation 
of the last term in the master equation gives (see Appendix A) 

N{J,9)=Y^XAJ)e m[0+n{J)t] N n {J), N n (J) = f J~o/"? ] ■ ■ (59) 

A — ' J Zir lo + nil [J ) — ie 



12 



E. Weak dissipation limit: averaging over angles 



So far, we have kept / arbitrary. To study tunneling, however, we may impose the additional condition that 
/ = f(J)> an( i obtain a simpler equation by averaging the Fokker-Planck equation over the angle variable 6. This 
approximation has been discussed by Kramers |7|69[| as prevailing in the weak dissipation limit. Recall that then 
{H s , /} = 0, and that, for any phase space function 9(6, J), 



j>d6 {X,9} = <j> dX 



{ 89 


dX 


i) 


\ 8J 


o 9 J 





dJ J dJ 



— <p dO p9 



where we have used that for the classical trajectory d6 = 6dt — fldt, and that dX = dgX\jdd + djX\gdJ from where 
by imposing dX = we can deduce 8j6\x- 

Finally the Fokker-Planck equation (jig) becomes 



9 -L = ±{N 
8t dJ ^ 



df_ 
QdJ 



Vf 



where we introduced T> and JV as follows 



v = - Id8 P r = iJ2\Xn(J)\ 2 nin(J), 

J n 

N = jdep^ = nY^\Xn(J)\ 2 n 2 N n {j). 



(60) 

(61) 
(62) 



Now observe that from Eq. ( |58| ) and using that l/(s + ie) = PV(l/s) + ittS(s) we can write 

duj loj (u) 



ln (J) = PV 



--n£l{J) 1 [nSl{J)), 
Z7T to + nil (J) 2 



but since 1/(uj — nil (J)) — l/(u> + nfl (J)) = 2nf2/(w 2 — n 2 tt 2 (J)) the first term of j n above integrates to zero, so 
that only the imaginary term contributes, and we finally have 



V = ^\X n (J)\ 2 n 2 -y[nn] 



(63) 



A similar computation using expression (B9J) for N n leads to the final expression for J\f: 



(64) 



F. A rough estimate of the decay rate 

Eqs. (Q) and ( |60| ) are the basic equations for the rest of our analysis. The rest of the paper is devoted to the 
explicit computation of the T> and J\f functions in a field theoretical problem, and to solving the dynamical equations 
therefrom. However, we may already at this point make an educated guess about the relationship between the decay 
rate predicted by these equations, and the usual quantum estimates. 

The point is that, since these are after all equations similar to those discussed by Kramers M, we may obtain a rough 



estimate of the transition amplitude by just plugging in Kramers' result for the flux. This is made from Eq. (60) which 
we may write as a continuity equation dtf + djK = where the probability flux K may be directly identified from the 
equation. Then one looks for stationary solutions with positive flux Kq, which must satisfy (Af/Cl)djf + T>f = —Kq. 

From this equation one may estimate (imposing that the particle is in the potential well J Js f(J)dJ < 1, where J s is 
J at the separatrix, i.e. the top of the potential barrier) using that dE = fid J the following upper bound for Kq 
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K ~ e- B(E °\ B = J dE (V/Af) . 

At high temperature, we have v = jkT, Af = kTV, B = E/kT, and 

K ~ exp (—E s /kT) , 

where E a is the energy at the separatrix, as expected p9[ |. At low temperature, see Eq. (|55|), v = %\u)\^fj2. Now, 
because the sums which define both Af and T> are d ominated by frequencies of the order of the curvature of the 
potential around the unstable fixed point to ~ \J \ V" | (this being the timescale for the exponential approach to the 
unstable fixed point), then Af ~ hy / \V : !'\'D/2, and thus 

K Q ~ exp (-2^ s /?i^T^) . (65) 

We may compare this estimate with the usual WKB tunneling amplitude given by 

K{){tunnel) ~ exp 

(-ft -1 J dx \Z 2 V( x )^j ■ (66) 

If the integral is dominated by the peak at the unstable fixed point X s , then V ~ E s — V" (X — X s ) 2 /2 and the 
Euclidean trajectory may be parametrized as X ~ A s — y / 2E s /V" cos 9, < 9 < n, which gives / da; ^/^^ (x) ~ 

We arrive at the remarkable conclusion that tunneling is comparable to the noise effect from the environment. Of 
course, our coarse estimates are not reliable for an accurate comparison and we must proceed to a more quantitative 
account. 

IV. QUANTUM FIELDS AS OPEN SYSTEMS 

A. The model and system-environment split 

We wish to study vacuum decay at zero and finite temperatures for a 3+1 quantum massive scalar field $ with 
action 



(67) 



where the metric convention is rf- v — diag(— 1, 1, 1, 1) (/i, v = 0, 1, 2, 3). Although we keep fi explicit, we set c = 1. 
The mass M has units of (length) -1 , i> has units of My/h and g of M/y/K. 

As discussed in the introduction, we wish to focus on the behavior of the long wavelength modes of the field. Let 
us split the field $ = <fi + cp, where <f> represents the long wavelength modes and ip the short wavelength modes. To do 
that we introduce a length scale A -1 which will be suitably fixed for our problem. To define <p(x) we take a window 
function W(x' — x) centered at the point x with a width A -1 and covolute the field $ with it 

« X ) = f**W« -,).(✓). (68) 

then, of course, tp(x) — &(x) — <f>(x). In this way the field at each point has two contributions one corresponding 
roughly to scales larger or of order A -1 and the other to smaller scales. This has its correspondence in momentum 
space, we may define the Fourier transform of $ by 

*(*) = j (|^3 e **(*)' ( 69 ) 

the long wavelength modes now become (f>tit) — W(— where W(k) is the Fourier transform of W, and for the 
short wavelength modes we have (ft = — 4>z- It may be convenient to use a Gaussian window in this way W is also 
Gaussian with a with A, or some times it may be more convenient to take a step function for the Fourier transform 
of the window such as W(k) — 9(A ~ k) where k = \k\. 
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At this point <f> is still a field, that is, it contains an infinite number of degrees of freedom. It is convenient to 



reduce the system to a single degree of freedom, such as we have discussed in section One way to accomplish this 
is simply to enclose the field in a box of size A -1 ; see 0. The boundary conditons in this case would introduce an 
undesired discrete spectrum for the modes, from a physical point of view it is more satisfactory to proceed as follows. 
In a region of volume A -3 we define the average field 

4>(t) = A 3 / d 3 x(t)(x). (70) 

Note that if we introduce the function p(k/A) — A 3 J A _ 3 exp(— ik ■ x), which satisfies p(0) = 1, then in momentum 
space we have <fi(t) = (2ir)~ 3 J d 3 kp(k/ A)0g so that 4> is made up of the modes of the field with k < A. 

Now when the the fields (j> and ip are substituted into Eq. (|S7|) the action will be decomposed in three parts. 
One involves the field (j) only, another the field ip only, both with the same functional dependence as the original 
action and a third interacting part involves terms linear in ip and the quadratic term (l/2)gipip 2 which comes from 
the cubic term in Eq. (|67|). We will approximate the different terms as follows: J A _ 3 d 3 x<f> n {x) ~ A 3 cj) n for any 
integer n and J A _ 3 d 3 xtp(x) ~ since the short wavelength modes should average to zero. Thus, generally we have 
Ja- 3 d 3 x(j) n f (ip) ~ 4> n J A - 3 d 3 xf{ip), where / is an arbitrary polynomic function of ip. Since we want to focus in the 
single degree of freedom <fi(t), which is the averaged field associated to a certain region of volume A~ 3 , it will be 
convenient to restrict the volume integrals involving the field <fi to that volume. 

Up to this point the scale A is arbitrary, now we need to fix it. Since the field 4> is homogeneous in region considered 
the gradient terms should be negligible in front of the mass term, (Vcf)) 2 <C M 2 (j> 2 and this means that the momentum 
cutoff for the modes of the field <j> should be k <C Al, which implies that there is an upper bound for the scale A: 
A < M. We also want to be able to treat the field <p perturbatively and that means that the field should be stable, in 
this case the cubic term ip 3 will be a two-loop order term and may be neglected to leading order in h. We will see in a 
moment that the condition for tp to be stable implies a lower bound for A: A > Al . Thus for the system-environment 
split to be consistent we need A ~ M and the volume is AT~ 3 . 

Finally the action ( |67| ) can be aproximated as a system-environment interaction action S[<f>] ~ S s [4>] + S e [ip] + 
Sint [0, <p] where 

SM = M- 3 [ dt (\$ )-M 2 4> 2 + Ig^ 3 ) , (71) 



2" 2 r 6* 



S e [(p] = dt 



A>iU 



d 3 k 
(2^ 



- ~ (fc 2 + m 2 ) ipj-cp^ 



(72) 



S int [<p,(p] = --g / dt [ d 3 x<j)ip 2 , (73) 
1 J Jm- 3 

where we made a mode decomposition and integrated over the whole space volume in the environment action. Now 
the potential of the system 

V$) = \m 2 4> 2 - i ff 3 , (74) 
2 o 

has a stable fixed point at <fi — and an unstable fixed point at 4> = 4> s = 2g~ 1 M 2 . The former corresponds to 
zero energy, and the later to E — E s = Af _1 0j/6, in the volume M~ 3 . For intermediate energies, we may have 
bound and unbound states. They are separated by a potential barrier, which at zero energy extends from = 
to 6 

— <&exit = 30 s /2 . At any given energy there will be three classical turning points cf>i, < < cf>ji < cf> s < <\>x\ 
as E — > 0, 4>l, 4>R ~^ an d 4*x —> i'exit, while when E — > E s , <j)R,4>x <Ps and <j)L — * —(f> s /2. We are thus in 



the situation described in Sec. ^I|, consequently according to the estimate at the end of last subsection, Eq. (|6£ 
we expect here that the tunneling rate will be Kq ~ exp [—aM 2 /Tig 2 ^ where a is a dimensionless parameter to be 
determined. Since the system is the averaged field cj> the tunneling rate is now per unit volume. When comparing 



with Sec. Ill note that the system coordinate x and frequency S7o in the QBM model correspond here to the averaged 
field (j) and the field mass M, respectively; also here there is only one coupling parameter A = g. 

For the environment modes to be stable as required we should have from the previous actions S e and Si nt that 
k 2 + Al 2 — g<f) > 0, and since the maximum value that <p may take at the barrier is <fi s = 2g~ 1 M 2 , it is clear that 
k > M and this gives the lower bound for A mentioned above. 



Let us now follow Sec. Ill and define the function S(t) which appears in the interaction action which is a key 
function to construct the dissipation and noise kernels, see S_y(t) in Eqs. ( p"3] ) and (|l4|), or E(t) in Eqs. (51) and (p32]). 
Thus we likewise define 
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S(t) = M 3 / d 3 4s<P 2 (x) , 



(75) 



and the interaction action becomes Smt = M 3 J dt<j)(t)E(t) ; note that S(i) has units of M 3 \fh. By analogy with 
( |l0| ) the influence functional becomes 



e iS IF /h _ I D(pDlp ' pe . exp . 



S e [p] ~ S e [<p'] - M~ 3 / di (0(t)3(i) - 0'(t)S'(t)) 



where is the environment density matrix at the initial time, 
quadratic terms and we can write 



(76) 



As in the previous section, see Eq. (jl5|) , we keep only 



S 



1 



IF 



~-M~ 



b')-H 



V) + - 2 M- 3 { 



■4t).N-$-4f). 



(77) 



Thus we are now in the situation described in general terms in Sec. Ill, and must complete the following steps: (a) 
identify the kernels 7 and v 1 (b) evaluate the Fourier transforms of x and p along a classical trajectory, (c) compute 
the functions V and AT, (d) solve the Fokker-Planck equation for /, or at least justify the approximations already 
discussed, and (e) evaluate the decay rate. We shall carry tasks (a), (b) and (c) in this section, and leave (d) and (e) 
for the next. 

We should emphasize that as remarked in the previous section, the kernel H needs to be regularized and that this 
involves a mass renormalization, we assume from now on that the mass M which appears in the action as well as 
the coupling parameter g are the renormalized values. Also we emphasize that we have chosen to reduce the low 
frequency part of the field to a single degree of freedom only for simplicity. An inhomogeneous field could be handled, 
for example, with the techniques presented in Ref. [[75[. 



B. The dissipation and noise kernels: 7 and v 

In order to find the noise and dissipation kernels, we use the representations in Eqs. (|5l]) and (^2|). To compute 
the averages, observe that the environmental variables can be treated as a free field. The details for the evaluation of 



the correlation (S(i)S(i')) are reproduced in Appendix B, see Eq. (109), where we find 



,„ , v „ , ,„ ,„ , v, ,„ , /\ 1 h 2 g 2 M 3 f d 3 p e -Mwp(*-*') 

(5 (t) 5 (0) = (S (t)) (S (0) + / —^3 -2 • (78) 

8 Jp>M (27Tj w p 

One may compute the anticommutator and commutator operator required for the kernels ( |5l| ) and (|5^ from this 
expression. Then comparing with Eqs. (|53|) and (|54]), and recalling that the dissipation kernel must be corrected by 
a factor of M -3 , in agreement with Eq. (|77|), we obtain 

ittfg 2 f d 3 p 



huj-f (w) = — — / - 3 [5 (2uj p ~u)-5 (2uj p + u)] 
4 J P >M (2tt) cj£ 



nh 2 g 2 . , . f°° 4irp 2 dp r . 
2— ^sign (u) / ^ 3 \ S (2u; p - \oj\) 

Ion V oj z 



which leads to 
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With a similar computation, for the noise kernel, we get 



W 4M 2 , 9 9 . 



, N h 2 g 2 M 3 I 4M 2 , „ , s 

" (w) = "W" V 1 " ^ (w 8M ) ■ (80) 

These kernels are related as required by the zero temperature fluctuation-dissipation theorem, see Eqs. (|55"|). We 
observe that these are the same as the kernels found in Ref. 13 in a different context. The analysis there substantiates 



our claim that particle creation, and backreaction thereof, are the main mechanism for dissipation and noise in this 
model. 
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C. Classical trajectories 



Let us recall the basic definitions. The form of the classical potential energy density Eq. (J7J) suggests writing 
4> = <fi s x — 2g~ 1 M 2 x so that V(<j>) — 4g~ 2 M 6 (x 2 jl — x 3 /3). Let us introduce also a dimensionless time r = tM and 
the classical action density for the system, which is defined as S s = M 3 S S , reads 



S s 



4M 5 



(It 



1 SdxY 

2 \dr) 



v{x) 



(81) 



where v(x) — x 2 /2 — x 3 /3. This action density includes dimensionless variables only the critical point is now x = 1, 
corresponding to dimensionless energy e s = 1/6. Physical energy densities are of course E = (4M 6 /<7 2 ) e = 6E s e. 




FIG. 1. Plot of the dimensionless potential v(x), as a function of the dimensionless variable x — (f>/(f) s . Note that it exhibits 
a metastable minimum at x — and an unstable maximum at x = x s , where the value of the potential coincides with the 
energy e s of the separatrix trajectories. The three turning points xl, xr and x x corresponding to any energy e < e s have also 
been represented. 

For any energy e < e s , we shall have three turning points xl < < xr < 1 < x x , see Fig. [l]. The classical orbits 
are expressed in terms of Jacobi elliptic functions with modulus n given by 

k 2 = XR ~ XL j (82) 

Xx XL 

which satisfies that k — > when e — > and K — > 1 when e — > e s . The actual expression, computed in Appendix C, is 

x = xl -\- (xr — xl) sn 2 Ar, (83) 
where A = \J (x x ~ xl)/6- The function sn has period AK [k], and sn 2 has half of that, so the dimensionless frequency 

is n = \tt/k. 

The Fourier coefficients of x(r), for are computed in Appendix D. They are given by 

3 ^ 2 (84) 



sinh (nCLK'/X) 



D. The dimensionless Fokker-Planck equation 



Let us first write the noise and dissipation kernels in terms of dimensionless variables. For the system Hamiltonian 
density, i.e. H s = M 3 H S , besides x_ and r (4> = 2g~ 1 M 2 x , t = t/M) we write also p = 2g^ 1 M 3 p, V = Ag^ 2 M 6 v(x) 
(V'((p) — 2g~ 1 M 4 v'(x)), thereby H s = (4g _2 M 6 ) h s . Note that this transformation is not canonical but preserves 
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the equations of motion, since pd(f> — H s dt = (4g~ 2 M 5 ) (pdx — h s dr), see Ref. j76|. The action variable also becomes 
J = (4<7~ 2 Af 5 ) j, while the angles are, of course, unchanged. Since the transformation is not canonical, it does not 
preserve phase volume: dpdcf) = Ag~ 2 M 5 dpdx. Therefore, we ought to define a new distribution function / (<f>,p, t) = 
( ff 2 /4M 5 )/(x,p,r). Then we find df/dt -> ( ff 2 /4Af 4 ) df/dr, and {B s ,f} -> ( 5 2 /4M 4 ) {h s , /}._. . Therefore, 
instead of Eq. (Eq), we have 



d 



2M 4 dp 



r(&j>,t)/- 



AM 5 



{ N >f} s , 



(85) 



Now let us recall the definitions of 7„ and N n in Eqs. ( |58| ) and (|5^) and let us introduce the new functions u> = Mu>, 
7 (w) = (hg 2 /M) 7 (a)) , 7„ (J) = % 2 7„ (j) and <^„ (J) = 2g- 1 M 2 x n {j) to get T = 2%M 2 f r). On the other 
hand, let us call v(u) = h 2 g 2 M 3 D(ui), N n = h 2 g 2 M 3 N n and N = 2h 2 gM 5 N and we finally obtain the following 
dimensionless form of the Fokker-Planck equation 



df 
dr 



4 d 

(3 dp 



(86) 



where (3 = AM 2 /(Tig 2 ). In dimensionless variables, the energy scale for the false vacuum is 1/(2/3). 

We may now proceed to write down the angle averaged Fokker-Planck equation ( |60| ) in dimensionless variables. 
This takes the form 



df 
8t 



4 d ( N df 



5/ 



(87) 



where D = (0/2) J2 n \x n (j)\ 2 n 2 j [nO] and N = (0/2) J2 n \ x n (j)\ 2 n 2 v [nfi] , and we have also introduced the 
dimensionless = 0/M. Using the explicit form of the Fourier coefficients Eq. ( p4[ ) and the expressions @ and 
( |so| ) for the dissipation and noise kernels, 7 and 1/, we have explicitly 



D 



N 



90 

16tt ^ sinh 2 (nfiif'/A) 

v ""' 

32tt ^ n sinh 2 

n>0 



n 2 2 



-6»(n^0 



2o2 



iflK'/X) 



n 2 2 



-6»(n^0 



2o2 



(89) 



There are two relevant limiting cases. The first limit corresponds to the separatrix energy, that is, when e — > e s , 
— > 0, A — > 1/2, if' — * 7r/2. We may write £ = nO and ^ n ~ (l/O) J d£ in the previous equations. Then after 
numerical integration we obtain 



D = D S 



16tt 



-2.441 ... x 10" 



N = N. 



9 

32tt' 



-7.381... x 10 



-7 



(90) 



which remain finite. The second limit corresponds to the bottom of the potential, that is, when e — > 0, O — > 1, A — > 
1/2. Then K' ~ (1/2) ln(16K- 2 ) = (1/2) ln[24(x/j - xl)' 1 ] ~ (l/4)_ln(72/e_), see Ref. Q, formula 8.113.3, and 
sinh -2 (nflK'/X) ~ 4 exp(— 4n_ft' / ) ~ 4(e/72)™ and the expressions for D and N take the values 



47T ^ 

n=3 



(4)'" 



6e 3 , N , 



E- 4 

n — 3 



(91) 



where a and b are numerical coefficients that can be read from these expressions. Note that the rapid decay of these 
functions as e — ► is due to the presence of a threshold, encoded in the theta functions of Eqs. ([79]) and (|80|), which 
enforce that n 2 2 > 8. At finite temperature, where the self-energy of the fluctuations remains complex even on-shell 
(when higher loops are included, of course), D and N fall like e, as in Kramers' original analysis M. 



V. ENVIRONMENT INDUCED DECAY RATES 



We may now conclude our task of finding the decay rates as described by the Fokker-Planck equation. In agreement 
with our weak dissipation assumptions we shall assume that we may average over angles, and restrict our analysis to 
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the Fokker-Planck equation (|S7|); a related analysis of angle-dependent solutions is given in Appendix 6 of Ref. J63| . 
Since the equation is linear in the Wigner function, it is convenient to first seek the normal modes, namely, solutions 
with a simple (exponential) dependence on time. The desired solution will then be reconstructed as a superposition 
of normal modes |39],([8) . 

Assuming then that / = e~ rT F (j), we get from (p7p the time independent equation 



LF + rF = 0, 



L-~— (— — 

~ 13 dj [f3(ldj 



D 



This equation may be written as a conservation equation LF = d<fr/dj, where $ is the flux defined by 

- 4 / N dF - 



S3 \pn dj 

Observe that for any e > e' > we have the following relationship between F, $ and r: 

$(£)-$(£') 



SI Fdj 



(92) 



(93) 



(94) 



Let us now introduce the unnormalized equilibrium solution Fb 
written as 



d 
dj 



n d 



F 
~Fb 



rF = 0, 



exp(-/3 / dj fiD/N), then Eq. (||) can be 



(95) 



where p 2 (j) = 4NFe/(fj/? 2 ). Multiplying by F/Fb and integrating we see that r must be non-negative, as expected. 
From the mathematical point of view, this is an eigenvalue problem of the Sturm-Liouville type, and we may handle 
it in the usual way p9j. Let us begin by analyzing the solutions with r ^ 0, which are the decaying solutions. 



A. Decaying solutions 



Let us seek a solution Fk of Eq. ( |95| ) with nonzero r. Since later on we shall be interested in the long time behavior 
of solutions, we may focus on the range of small values of r. For concreteness, let us assume r[3 3 /3a -C 1. Let us 
introduce h by F K = {h/p)F B = {!3h/2)^j£lF B /N and write Eq. (§|) and the flux as 



d_ 

dj 



d (h 



9 d J \P 



rn(3 2 



ph = 0, 



9 d 
dj 



Expanding the second derivatives dj[p 2 dj(h/ p)] = dj(ph' — hp') = ph" 
respect to j, we obtain 



h 
P 

hp" 



{ph 1 -hp'). 



where dj stands for the derivative with 



h" 



rn/3 2 
4N 



)h = 0. 

P 



(96) 



We have two regimes in this equation: for large e (but not exponentially close to s s = 1/6), p is dominated by Fb, 
and p" j p ~ /3 2 fi 2 D 2 /(4N 2 ) dominates the r term, which is negligible. For e — * 0, on the other hand, p ~ e 3 ^ 2 , 
p" I p ~ (3/4) e~ 2 , and the r term is the leading term. The transition occurs for e — e* ~ r(3 2 /(3a), which by 
assumption is much smaller than the energy scale 1/2/3 typical of the false vacuum. 

For large e, the equation h" — (p" / p)h = admits of course the solution h = p, which gives back the equilibrium 
solution, i.e. Fk = Fb- To find the second solution we may write h — K r pa (here K r is a constant) to get 
a" j a' = —2p'/p which implies a' = —p~ 2 . This second solution corresponds to a constant flux $ = K r . and also is 



2 . This second solution corresponds to a constant flux $ = K r 
the negative exponential WKB solution. For small e, the equation h" + (r/3 2 /4ae 3 )/i = has solutions 



h = Cry/eZx 



(97) 



where k 2 = r(3 2 /a, Z is a Bessel function and C r is a constant (observe that all solutions are bounded). As e — > 0, 
F - e~ 3 / 4 cos (fce^ 1 / 2 ) but $ - e 3 / 4 sin (fe- 1 / 2 ) -> 0. 

Note that we obtain solutions for arbitrarily small values of r. This behavior, which is unlike that found by Kramers 
0, reflects the existence of a threshold for dissipation: at arbitrarily low energies, the motion is essentially harmonic, 
there is no dissipation, the fluctuations switch off (in agreement with the fluctuation-dissipation theorem) and the 
system requires an infinite time to climb out of the potential well. 
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B. Hilbert space structure and normalization 



The structure of Eq. (^5|) suggests introducing an inner product 

(f,9)=£^f*9, (98) 

where the star means complex conjugation. Then, if F r and F s are the solutions corresponding to eigenvalues r and 
s, we can write (r-s){F s ,F r ) = - f^ 3 dj F^ l [F s {LF r ) - [LF s )F r ] = -F^lF^r - $ 5 -F r ]|g S . Imposing the boundary 
conditions F r (j s ) — 0, we get rid of the contribution from the upper limit. In the lower limit, we may use the 
asymptotic form of the Bessel functions to get [F s $ r — $ s F r ](e) ~ 2n~ 1 C' r C' s k r 1 ^ 2 k^ 1 ^ 2 sin[(/c,. — k s )e~ 1 ^ 2 ] which 
converges (weakly) to zero. It is therefore natural to adopt the continuum normalization prescription 

(F s ,F r )=S{r-s). (99) 

As in Landau and Lifschitz' analysis of the WKB wave functions [||, this singular behavior is caused by the 
oscillations as e — > 0. More precisely J Q dj F I3 1 F r F s ~ ir~ 1 C r C s k r ~ 1 ^ 2 kJ 1 ^ 2 L de £ _3 / 2 cos[(fc r — k s )e~ 1 ^ 2 ] ~ 
4C 2 /c J T 1 <5(fc r — k s ) ~ 8aC 2 f3^ 2 5(r — s), where the integral upper limit is anything. We thus find the constants 

C r ~ /3/VSa. (100) 

To find the constants K r we should match our solutions accross the transition point at e*. Without getting into 
details, see Ref. p0| |, it is clear that in order of magnitude 

i nr/fe^ 3 3 / 4 

r ~ p(e*)tr(e*)V irk r ^ ~ 40Fr x /V ( E *) £* 3 / 4 ' ( ' 

To find cr(e*), we recall that according to our boundary conditions cr(j s ) — 0. Therefore c(e*) = j e ° djp~ 2 (j). 
Let us split the range of integration in the two segments [e*, 1/2/3] and [1/2/3, e]. The contribution Ijj from the 
second segment is dominated by the exponential growth of F^ 1 near the separatrix, thus In ~ FB(l/2/3)Fg 1 (e s ) ~ 
A' -1 exp(/3D s /6N s ) where A' ~ O(l). In the first segment, we find p~ 2 — /3 2 /(4ae 3 ), so it contributes as Ij ~ 
(l/8)/3 2 a- 1 [(£*) -2 - (V4)/?- 2 ] ~ (l/8)/3 2 a- 1 [9a 2 /3" 4 r- 2 - (l/4)/3~ 2 ] ~ (9/8)af3- 2 [ r - 2 - (l/36)/3 2 a^ 2 ]. Adding both 
contributions, and keeping only the leading terms, we find 

I *\ 9a / 1 1 
where A ~ A exp(— /3D S /12N S ) with A ~ O (1) . So finally we have 

K r~K^#> (102) 

where K is a r-independent constant. 



C. Decay of the false vacuum 



After all this work, the dynamical problem is now trivial. We are interested in a time-dependent solution with an 
initial condition similar to the Wigner function of the false vacuum, F = Fpv ~ 2/3 exp (— 2(3e) 0|. The solution is 



poo 

F(j,t)= / dr 
Jo 



-rt 



c r F r (j) 



where c r = Jq S dj Fg 1 F r Fpy . Fortunately, we are interested in the range of small r, where F r peaks at values much 
smaller than 1/2/3. The only feature of Fpy and Fb that we really need is that they are smooth there. Thus, 

c r ~ FfvWFb 1 ^) ft $ F r = WKr/r, (cfr. Eq. Q). Finally 
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/■oo -rt 

F (j, t) = 2(3K\ 2 J dr ^~^Fr (j) ■ 
To find the persistence probability, P(t), we integrate F(j, t) over the potential well. Thus after using the previous 



result for J g 3 dj F r (j) and the value of K r from (102) we obtain 

AV 



P (t) ~ 2(3K 2 / dr 
Jo 



fr 2 + A 2 ! 2 ' 



The analytic expression for this integral is given in Ref. 0, formula 3.355.1. When Xt ~ 1 it can be approximated 
with an integral, the best fit is obtained when P (t) ~ exp(— 0.4At). For larger times, we have a crossover from 
exponential to power law 1/t decay, as expected from quantum mechanics (cfr. p5| ). To summarize, we have proven 
that the Fokker-Planck equation leads to an exponential decay of the false vacuum, with a decay rate t^ 1 ~ A: 

^ 1 -Acxp( I ^ /3 ), (103) 

with A ~ 0(1). This equation may be compared to Eq. (fi| ) , here (3 — AM 2 /(kg 2 ), as defined in Eq. (|86|), and the 
ratio T> s /N s = D s /N s where D s and N s are given in (^)l[); see Eqs. (p^), ( |64| ) and (|^) for the definitions of the 
functions involved. 



D. Comparing with the instanton method 



Let us conclude by comparing the rate estimate of Eq. (103) with the estimate derived from the instanton method. 
The Euclidean action for our model is 



S E [$] = j d A x fy^dpQdv* + \m 2 ^ 2 - l -g^ , 



(104) 



where 5^ stands for the Euclidean metric. We are interested in SO (4) symmetric instantons, which depend on all 
four Euclidean coordinates but only through the Euclidean radius p = {x 2 +y 2 + z 2 +t 2 E ) 1 / 2 , where tg is the Euclidean 
time. Let us scale $ = <fi s f (Mp) to get the Euclidean action (in d dimensions) Sd = (4tt/3) h(3Id, where 



Id = dp p' 



„d-l 



(105) 



Of course, computing the four-dimensional (d = 4) instanton is not simple, but we may approximate I4 by ii, in 
which case we may use the simpler one-dimensional formula 



S 



3/2 



dx \/2v(x) 



dx \ x 2 
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All in all, the instanton prediction is lnt rf 1 ~ —1.2 . . . x (3. On the other hand, the noise induced prediction, Eq. ( 103 ), 
is hii^ 1 ~ In A ~ (— D s /12N s )/3 ~ —0.05 . . . x (3. So in this simple case the noise induced contribution dominates 
over the instanton contribution. 



VI. CONCLUSIONS 



In this paper we have studied the contribution to vacuum decay in field theory as a consequence of the interaction 
between the long and short-wavelength modes. We have seen that the dynamics of the long-wavelength modes 
becomes diffussive in its interaction with the short wavelength modes. On the one hand, there is dissipation of the 
long-wavelength modes due to the excitation of the short-wavelength sector, and in turn, that latter sector induces 
fluctuations into the first sector. As a result there is a significant contribution to the total decay rate due to activation, 
even at zero temperature. 

A few remarks are in order. What we have shown is that if we consider the field initially in a mctacstblc phase 
in a region of size M^ 1 there is a probabiliy per unit volume and unit time of decay to a stable phase given by 
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Eq. (103). Of course, the size of the bubble formed is of order M _1 , in our calculation this size is fixed. We cannot 
consider smaller size regions because we could have not neglected the gradient terms of the system in front of the mass 
terms, also we could not have considered larger size regions because the environment modes would become unstable, 
their evolution would become nonlinear and our perturbative treatment of the influence functional would break down. 
However, the critical bubble size that one obtains in first order phase transition in statistical physics or in field theory 
using instanton methods is also of the order of M . This, in our opinion, makes this computation of interest, since 
one expects that once the critical bubble is formed it will evolve in the usual way; see for instance Ref. ||. Had our 
calculation involved a size smaller than the critical size then the bubbles formed could not grow and would collapse, 
as then the energy of the bubble wall would overcome the energy difference between the metastable phase and the 
bubble interior. 

Another relevant point concerns the energy balance in the activation mechanism of vacuum decay described here 
which involves noise and dissipation. In Appendix E we show that the average power exchanged between system and 
environment is zero: if some trajectories gain power through noise, some others lose power through dissipation. Of 
course, the balance is only statistical, but we must stress that this gain and loss process would go on even if there 
were no separatrix and the system were in equilibrium. The only reason why the system does not equilibrate in our 
problem is that we remove those particles that reach the separatrix, as demanded by our boundary condition there. 
We should say that the system equilibrates but we assume that the true vacuum is very much deeper than the false 
vacuum, so the equilibrium distribution vanishes inside the potential well. Note that this is an expected result in 
vacuum bubble formation. Once the critical bubble is formed the energy released in the conversion from false to true 
vacuum is converted into energy of the growing bubble wall, so that the energy balance is still zero, this last aspect 
however cannot be studied with the present analysis. 
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VIII. APPENDICES 



A. Appendix A 

Our problem is to compute an expression like 

j dt'N(t-t')(^^—R[X(t),p(t)}\, (106) 

where R is an arbitrary functional of £(i) and (...) denotes expectation value with respect to Pq[£;t). Here, for 
simplicity, we will assume also that Pq is independent of X as in the case of the cubic potential. In our case, this 
leads to 

4 H . N <« - (f$« p ) + 1 . 

where S x = 6{X(t) - X) and 8 P = 6{p(t) - p). 

Within the "reduction of order" procedure, we substitute dX (r) jdX it) and dX (r) /dp (t) for SX (t) /S£ (t) and 
Sp (t) /<5£ (t) , where the variations are understood in the sense of the result of coupling a stochastic source to the 
classical equations of motion. The argument runs as follows: we know that SX (r + ) /<5£ (r) = and Sp (r + ) (t) = 1. 
On the other hand S ((t) X(t+) = d x{t) X( T +)S s{T) X(t) + d p{t) X(r+)S c(T) p(t) and % r) p(r+) = d x{t)P (r+)S s{T) X(t) + 
9 P (t)P( T+ )^(T)P(t)- By Liouville's theorem the determinant of this two by two system is one, so indeed 
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SX(t) _ dX{r+) 6 P (t) _ 8X(t+) 
SUr) ~ dp(t) ' 6^) ~ dX(t) ■ 
The right-hand side of these equations are continuous, so we may omit the superscript. Furthermore 



dX(r) 



= {X(t),X(r)}, 



dX(r) 



= -{p(t),X(r)} 



dp(t) w ' WJ ' dX(t) 

and computing the Poisson brackets in terms of the canonical variables 9 and J we arrive at 



= h dt'N(t-t') 



A [{X (i) , X (f )} /] + ^ [{p (t) , x (*')} /] 

= ^ [{* (*) , ^ ( J, 61)} /] + ^ [{p (t) , N (J, 0)} /] , 

where we have defined N (J, 9) = Ti J dt' N (t — t') X (t 1 ). Applying the reduction of order procedure by substituting 
X(t') by the Fourier expression (p6|) and using Eq. (p4|) the function N(J,0) is written as, 



N (J, 0) = X n (J) e m[e+n{J)t] N n (J) . 



(107) 



where N n (J) is given by Eq. (|59|). We may simplify further, by observing that {X, {p, N}} — {p, {X, N}} = 
— {N, {X,p}} — by the Jacobi identity and the fact that {X 7 p} — 1, then by further manipulation of the Poisson 
bracket terms we finally get 



m = {X, {p, N (J, 9)} /} - {p, {X, N (J, 9)} /} 
= {p, N (J, 9)} {X, /} - {X, N (J, 9)} {p, /} 
dpdX dXdp~\ \dNdf dN df 
dJ~d9 + ~dJd9 ~d9dJ ~ ~dJ "d9 



= -{NJ} 



(108) 



B. Appendix B 

We first expand the environmental variables in creation and annihilation operators 



h 

2u> k 



a k e 



where lj 2 = k 2 + M 2 . Next, recalling the definition of 3 given in Eq. rt75|) we may write 



Z(t) = -gM 3 



M-3 



d 3 k 

(2^ 



d 3 p 



P ,k-p>M (2ir) 



To perform the integral over x we introduce as in Sec. IV the function p (k/M) — M 3 J M _ 3 d 3 x exp(ifc • x), which 
satisfies p(0) = 1. Then by direct substitution we can write the correlation (Z(t)Z(t')) as 
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(3 (t) 5(f)) =(S(t))<H(f)) + 



d 3 k d 3 k' 
(2tt) 3 (2nf 



I- 

P\ 77 \P\ 77 



M 



M 



f J^L f 

Jp,k-n>M (2tt) 3 Jp> 



d 3 p' e ~ iuJpt e i ^p'- k ' t ' e~ ibJk ^ pt e iiAj p' t ' 



a p a k - p al p ,al,_ k , 



lp,k-p>M (2tt) Jp',k'-p'>M (2n) y/UJpUlp-kUJp'LOp' -k' 

Now the vacuum expectation value in the last equation can be written as 

(apa k - p al p ,al,_ k ) = (2tt) 6 6^ (fc + fc') {<5< 3 > (k-p + f) +6^ (?+?)}, 
and since in any case we deal with values of k and kl much lower than typical values of p and p' , we get 



(3(t)S(f)) = <S(i))(S(f)) 



h 2 g 2 f d 3 k 



d 3 P e -2M*-0 

MJJ p > m (2tt) 3 u>1 



(109) 



To obtain Eq. @, observe that / d 3 k p 2 (k/M) = (2n) 3 M 3 . 
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C. Appendix C 



Let us write the integrand in the action density ( |Bl| ) as 

e - v(x) = -(x- x L ) (x - x R ) (x-x x ). (110) 

It can be rearranged as e — v(x) = (l/3)(x x — x)[(1/4)(xr — xl) 2 — (x — (1/2)(xr + xl)) 2 ], which suggests writting 
x = (1/2)(xr + xl) — (l/2)(a;fl — xl) cos(2</j) — xl + (xr — Xl) sin 2 ip, and thus 

e - v(x) = — (x x - x L ) (x R - x L ) 2 sin 2 2(p (l - k 2 sin 2 ip) , (111) 

where k is defined in Eq. (|S2|). 

The equation for a classical trajectory is r = J dx p^ 1 = J dx [2(e — v(x))]^ 1 ^ 2 and since dx = (xr — xl) sm(2ip)d(p 
we have 



6 r v (112) 



J o yl — k 2 sin 2 



When goes from to n/2, the point describes a half orbit. So the period T is twice the result from (112) when 
ip = 7r/2, and one can check that when e — * 0, T — * 27r, as it should. 

Using now the identities 16.1.3, 4 and 5 from Ref. plj] , we get the desired formula (|S3"|). Observe that as e — * 0, 
sn u ~ sinw, A ~ 1/2, and x = + (x^j — xl) sn 2 (r/2) ~ (1/2) (xr + xl) — (1/2) (xr — Xl) cost which corresponds 
to the harmonic trajectory as expected near the bottom of the potential. 

D. Appendix D 

Here we shall borrow an argument from Whittaker and Watson [Q. The Fourier coefficients of the classical 
trajectory are 

z n = i £ dre- inflr i(r) 



= x L 5 n0 + (x R - x L ) [ du e~ 2m ™sn 2 2Ku, 
Jo 



where we have changed to the u variable, used the explicit trajectories ( |83| ) and again redefined the integration variable 
in the second line. According to 16.7 and 16.8 in Ref. (83), we have the following properties: sn 2K (u + 1) = — sn 2Ku, 
sn 2K (u + 2iQ) = sn 2Ku, where Q = K'/(2K) with K' = K[yl — k 2 ]; and sn 2Ku has a simple pole at u = iQ, 
with residue 1/(2kK). 

To compute the Fourier coefficient (for n / 0) we consider the integral over the anti-clockwise contour T with 
vertices at 0, 1, 2iQ and 2iQ — 1. The two oblique sides cancel each other, and 



/ dz e- 2m * z sn 2 2Kz = q 4n / ds e- 2m7TZ sn 2 2Ks 

J-l+2iQ J -I 



fl 

= q- 2n I du e~ 2m ™sn 2 2K«, 



where q = exp(—irK'/K). Near the pole, we have (see 16.8 and 16.3.1 in Ref. |p3[| ) sn 2K (u + iQ) = sn (2Ku + iK') — 
(resn 2Ku)~ 1 = (2kKu)^ 1 [l + O (u 2 )] , and from Cauchy's theorem the integral over the contour T becomes 
2iri(4K 2 K 2 y 1 d z exp(-2imrz)\ z=i0 = wk 2 kT 2 K~ 2 q~ n which yields Eq. (||). 
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E. Appendix E 



Here we compute the average power exchanged between the system and the environment. Let us go back to the 
Langevin equations ( f49| ) with the Gaussian source described by Eq. (50). From these equations it follows that at 
the phase space point (X, P) there is a dissipated power Wd = —LP, where r = — f_ dt'H(t — t')X(t') and a noise 
power w n = £P. The power dissipated in the whole ensemble is Wd = — J dXdP f (X, P) TP. 

Let us use the action-angle variables at t = as Lagrangian coordinates identifying a given trajectory. Then the 
power dissipated is 



W d 



dJ f (J) / dO TP 



dJ f(J)n(J)V. 



The average noise power in the whole ensemble is (W n ) = J Q B dJ f (J) J d6 (£ (t) P) where P(6,J,t) has been 
perturbed away from the classical value by the action of the noise. We then get, using the Novikov trick 



dtN(t - t) WW) 



dtN{t ~ t) ^m 



{P,N}. 



We can simplify this by using the same arguments as in Appendix A: J dd{P,N} — J dd(deP)j(djN)p = 
dj J d9(d e P) ; N. Performing now a last integration by parts / dd {P, N} = dj J d9 N(d s P)j = -dj f dQ P(d 9 N)j = 
-djN. 



The total average power is thus 

W = W d + W n = 



dJf(J) 

Integrating by parts, knowing that TV" (0) = / (J s ) — 0, we get W — J^ s dJ D, (J) $ (J) where $ is the flux. Recalling 

now that £1 = dE/dJ, with our boundary conditions we have W = E S $(E S ) + J^ s dJ Edtf. We can now show 
that W ~ (this is obviously true for the equilibrium solution Fb, when both terms in the expression for W 
vanish independently). Let us write the general solution for the distribution / as the following superposition / = 



Since the flux is linear in /, and the flux for F r is K r , we get for W, 

rJ, 



w 



dr e rt c r 



E s K r 



dJ EF r 



The integral in the second term is dominated by the upper limit (since a ~ E 2 , the integral depends only logarith- 



mically on the peak E* ) , and L s dJ EF r 
as expected. 



E s J Q " dJ F r = E s K r /r, which makes the total averaged power W ~ 



[1] A. Guth, in The Very Early Universe edited by G. Gibbons, S. W. Hawking and S. T. C. Siklos (Cambridge University 

Press, Cambridge, England, 1983), p. 171. 
[2] E. Shuryak, Phys. Rep. 61, 71 (1980). 

[3] P. Shukla, A. Mohanty, S. Gupta and M. Gleiser, Phys. Rev. C 62, 054904 (2000). 

[4] G. Gamow, Z. Phys. 51, 204 (1928); R. Gurney and E. Condon, Phys Rev 33, 127 (1929); N. S. Krylov and V. A. Fock, 

Zh. Eksp. Teor. Fiz. 17, 93 (1947); V. A. Fock, Zh. Eksp. Teor. Fiz. 42, 1135 (1962) [Sov. Phys. JETP 15, 784 (1962)]. 
[5] L. Landau and E. Lifshitz, Quantum Mechanics (Non Relativistic Theory) (Pergamon, London, 1958). 
[6] A. Migdal, Qualitative Methods in Quantum Theory (Addison- Wesley, New York, 1977). 
[7] H. Kramers, Physica VII, 284 (1940). 
[8] J. Langer, Ann. Phys. 54, 258 (1969). 

[9] S. Coleman, Phys. Rev. D 15, 2929 (1977); C. Callan and S. Coleman, Phys. Rev. D 16, 1762 (1977); S. Coleman, V. 
Glaser and A. Martin, Comm. Math. Phys 58, 211 (1978); S. Coleman and F. De Luccia, Phys. Rev. D 21, 3305 (1980); 
S. Coleman, Aspects of Symmetry (Cambridge University Press, Cambridge, England, 1985). 



25 



P. Hanggi, P. Talkner, and M. Borkovec, Rev. Mod. Phys. 62, 251 (1990). 

I. Affleck, Phys. Rev. Lett 46, 388 (1981); A.D. Linde, Nucl. Phys. B 216, 421 (1983). 

J. Garriga, Phys. Rev. D 49, 5497 (1994). 

E. Calzetta and B. L. Hu, Phys. Rev. D 49, 6636 (1994). 

J. Gunton, M. San Miguel and P. Sahni, in Phase Transitions and Critical Phenomena, Vol. 8 edited by C. Domb and J. 
Lebowitz (Academic Press, London, 1983), p. 267. 

M. Castagnino and R. Laura, Phys. Rev. A 56, 108 (1997); D. Arbo, M. Castagnino, F. Gaioli and S. Iguri, Physica A 
277, 469 (2000). 

J. Langer, Ann. Phys. 41, 108 (1967). 

N. G. van Kampen, Phys. Rev. 135, A362 (1964); J. Lebowitz and O. Penrose, J. Math. Phys. 7, 98 (1966). 
H. Suzuki and H. Yasuta, Phys. Rev. D 57, 2500 (1998). 

D. Bes and J. Kurchan, The Treatment of Collective Coordinates in Many-Body Systems (World Scientific, Singapore, 
1990). 

T. Banks, C. Bender, and T. Wu, Phys. Rev D 8, 3346, 3365 (1973). 
M. Gleiser, G. Marques and R. Ramos, Phys. Rev. D 48, 1571 (1993). 

A. Goncharov and A. Linde, Fiz. Elem. Chastits. At. Yadra 17, 837 (1986) [Sov. J. Part. Nucl. 17, 369 (1986)]; A. Linde, 
in Three Hundred Years of Gravitation edited by S. Hawking and W. Israel (Cambridge University Press, Cambridge 
(England), 1986). 

V. Rubakov, Nucl. Phys. B 245, 481 (1984). 

T. Tanaka, M. Sasaki and K. Yamamoto, Phys. Rev. D 49, 1039 (1994); K. Yamamoto, T. Tanaka and M. Sasaki, Phys. 
Rev. D 51, 2968 (1995); T. Hamazaki, M. Sasaki, T. Tanaka and K. Yamamoto, Phys. Rev. D 53, 2045 (1996). 
L. Mersini, Phys. Rev. D 59, 123521 (1999). 
A. Schmid, Ann. Phys. 170, 333 (1986). 



G. Bonini, A. Cohen, C. Rebbi and V. Rubakov, Phys. Rev. D 60, 076004 (1999); |yiant-ph/9901062 . 
J. P. Paz and S. Sinha, Phys. Rev. D 44, 1038 (1991); ibid 45, 2823 (1992). 

D. Boyanovsky and C. Aragao de Carvalho, Phys. Rev. D 48, 5850 (1993); D. Boyanovsky et. al., Nucl. Phys. B 441, 595, 
609 (1995). 

J-L. Gervais and B. Sakita, Phys. Rev. D 16, 3507 (1977); H. De Vega, J-L. Gervais and B. Sakita, Nucl. Phys. B 139, 20 
(1978); ibid 142, 125 (1978). 

I. Affleck and F. De Luccia, Phys. Rev. D 20, 3168 (1979). 
U. Weiss and W. Haeffner, Phys. Rev. D 27, 2916 (1983). 

W. Buchmuller, Z. Fodor, T. Helbig and D. Walliser, Ann. Phys. 234, 260 (1994). 
K. Schleich, Phys. Rev. D 32, 1889 (1985). 



G. Munster, A. Strumia and N. Tetradis, cond-mat/0002278 



V. Kiselev and K. Selivanov, Pis'ma Zh. Eksp. Teor. Fiz. 39, 72 (1984) [JETP Lett. 39, 85 (1984)]. 

S. Weinberg, The Quantum Theory of Fields II: Modern Applications (Cambridge University Press, Cambridge, England, 
1996). 

F. Wegner and A. Houghton, Phys. Rev. A 8, 401 (1973). 
A. Strumia, N. Tetradis and C. Wetterich, Phys. Lett. B 46 7 279 (1999). 



A. Strumia and N. Tetradis, Nucl. Phys. B 542, 719 (1999); |hep-ph/9904357| v2. 
M. D'Attanasio and M. Pietroni, Nucl. Phys. B 472, 711 (1996). 

S.-B. Liao, J. Polonyi and D. Xu, Phys. Rev. D 51, 748 (1995); S.-B. Liao and M. Strickland, ibid 52, 3653 (1995). 
D. Dalvit and F. Mazzitelli, Phys. Rev. D 54, 6338 (1996). 



K. Aoki, A. Horikoshi, M. Taniguchi and H. Terao, |iep-th/981205C 



E. Calzetta and B.-L. Hu, Phys. Rev. D 55, 3536 (1997). 
H.B. Callen and T.A. Welton, Phys. Rev. D 83, 34 (1951); W. Bernard and H.B. Callen, Rev. Mod. Phys. 31, 1017 (1959); 
M.S. Green, J. Chem. Phys. 19, 1036 (1951); R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957); Rep. Prog. Phys. 29, 255 (1966). 
A.O. Caldeira and A.J. Leggett, Ann. Phys. 149, 374 (1983); Phys. Rev. Lett. 46, 211 (1981); A.J. Leggett et al, Rev. 
Mod. Phys. 59, 1 (1987). 

U. Weiss, Quantum Dissipative Systems (World Scientific, Singapore, 1993). 
A. Bray and M. Moore, Phys. Rev. Lett. 49, 1545 (1982). 

K. Fujikawa, Phys. Rev. E 57, 5023 (1998); K. Fujikawa and H. Terashima, Phys. Rev. E 58, 7063 (1998); K. Fujikawa, S. 
Iso, M. Sasaki and H. Suzuki, Phy s. Rev. B 46, 102 95 (1992). 
S. Matsumoto and M. Yoshimura. ^ep-ph/0008025| . 



R. Bruinsma and P. Bak, Phys. Rev. Lett. 56, 420 (1986). 
W. Horsthemke and R. Lefever, Noise-induced transitions (Springer, Berlin, 1984). 
M. Griffoni and P. Hanggi, Phys. Rep. 304, 229 (1998). 

J. Zinn- Justin, Quantum Field Theory and Critical Phenomena (Clarendon Press, Oxford, 1989). 

D. Hochberg, C. Molina-Paris, J. Perez-Mercader and M. Visser, Phys. Rev. E 60, 6343 (1999); Physica A 280, 437 (2000); 
J. Stat. Phys. 99, 903 (2000); Phys. Lett. A 278, 177 (2001). 



2G 



M. Eguia and G. Mindlin, Phys. Rev E 61, 6490 (2000). 

A. Starobinsky, in Field Theory, Quantum Gravity and Strings edited by H. De Vega and N. Sanchez (Springer, Berlin, 
1986), p. 107. 

E. Calzetta and B.-L. Hu, Phys. Rev. D 52, 6770 (1995). 

C. Kiefer, D. Polarsky and A. Starobinsky, Int. J. Mod. Phys. D 7, 255 (1998). 
M. Mijic, Phys. Rev. D 42, 2469 (1990). 

B. -L. Hu and A. Matacz, Phys. Rev. D 51, 1577 (1995); B.-L. Hu and S. Sinha, ibid 51, 1587 (1995); A. Campos and 
E. Verdaguer, ibid 53, 1927 (1996); B.-L. Hu, Int. J. Theor. Phys. 38, 2987 (1999). R. Martin and E. Verdaguer, Int. J. 
Theor. Phys. 38, 3049 (1999); Phys. Lett. B 465, 113 (1999); Phys. Rev. D 60, 084008 (1999); ibid 61, 124024 (1999). 
E. Calzetta and E. Verdaguer, Phys. Rev. D 59 083513 (1999). 

A. Vilenkin, Phys. Lett. 117B, 25 (1982); Phys. Rev. D 27, 2848 (1983); 30, 509 (1984). 

R.P. Feynman and F.L. Vernon, Ann. Phys. (NY) 24, 118 (1963); R.P. Feynman and A.R. Hibbs, Quantum Mechanics 

and Path Integrals (McGraw-Hill, New York, 1965). 

A. Roura and E. Verdaguer, Phys. Rev. D 60, 107503 (1999 ). 



E. Calzetta, A. Roura and E. Verdaguer, quant-ph/0011097. 
B.L. Hu, J.P Paz, and Y. Zhang, Phys. Rev. D 47, 1576 (1993). 
H. Risken, The Fokker-Planck Equation (Springer-Verlag, Berlin, 1989). 

W.H. Zurek, Physics Today 44, 36 (1991); J.P. Paz, S. Habib and W.H. Zurek, Phys. Rev. D 47, 488 (1993). 
M. Gell-Mann and J.B. Hartle, Phys. Rev. D 47, 3345 (1993). 

E. Wigner, Phys Rev 40, 749 (1932). 

J.M. Sancho and M. San Miguel, Z. Phys. B 36, 357 (1979); 43 361 (1981). 
T. D. Lee and G. Wick, Phys. Rev. D 9, 2291 (1974). 

F. Liu and G. Mazenko, Phys. Rev. B 44, 9185 (1991); F. Liu and G. Maze nko, Phys. Rev. B 4 5, 6989 (1992); F. Liu and 



G. Mazenko, Phys. Rev. B 46, 5963 (1992); R. Wickham and G. Mazenko, |cond-mat/9609163| and 9607152v2; A.J. Bray, 
Adv. Phys. 43, 357 (1997). 

V. I. Arnold, Mathematical Methods of Classical Mechanics (Springer, Berlin, 1978). 

I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals Series, and Products (Academic, New York, 1979). 

C. W. Gardiner, Handbook of Stochastic Methods (Springer-Verlag, Berlin, 1994. 

R. Courant and D. Hilbert, Methods of Mathematical Physics (Interscience, New York, 1962). 

J. Heading, An Introduction to Phase- Integral Methods (Methuen and Co., London, 1962). 

S. Habib, Phys. Rev. D 42, 2566 (1990); M. Berry, in Chaotic Behavior of Deterministic Systems edited by G. Iooss, 
R. Helleman and R. Stora (North Holland, Amsterdam, 1983); A. Ozorio de Almeida, Hamiltonian Systems: Chaos and 
Quantization (Cambridge University Press, Cambridge, England, 1988). 
E.A. Novikov, Sov. Phys. JETP 20, 1290 (1965). 

Handbook of Mathematical Functions, edited by M. Abramowitz and I. Stegun (Dover, New York, 1978). 

E. Whittaker and G. Watson, A Course of Modern Analysis (Cambridge University Press, Cambridge, England, 1940). 



27 



